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Abstract 

We analyze the properties of the energy landscape of finite-size fully connected p- 
spin-like models whose high temperature phase is described, in the thermodynamic 
limit, by the schematic Mode Coupling Theory of super-cooled liquids. We show 
that finite-size fully connected p-spin-like models, where activated processes are 
possible, do exhibit properties typical of real super-cooled liquid when both are near 
the critical glass transition. Our results support the conclusion that fully-connected 
p-spin-like models are the natural statistical mechanical models for studying the 
glass transition in super-cooled liquids. 
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1 Introduction 



In recent years a significant effort has been devoted to the understanding of 
glass-forming systems. Recent theoretical and numerical results clearly show 
that the slowing down of the dynamics near the glass transition is strongly 
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connected to the potential energy landscape geometry. The trajectory of the 
representative point in the configuration space can be viewed as a path in 
a multidimensional potential energy surface [1]. The dynamics is therefore 
strongly influenced by the topography of the potential energy landscape: local 
minima, barriers heights, basins of attraction an other topological properties 
all influence the dynamics. 

The potential energy surface of a super-cooled liquid contains a large number 
of local minima, called inherent structures (IS) by Stillinger [2]. All states that 
under local energy minimization will flow into the same IS define the basin of 
the IS (valley). With this pictures in mind the time evolution of the system 
can be seen as the result of two different processes: thermal relaxation into 
basins (intra-basin motion) and thermally activated potential energy barrier 
crossing between different basins (inter-basin motion). When the temperature 
is lowered down to the order of the critical Mode Coupling Theory (MCT) 
temperature Tmct the inter-basin motion slows down and the relaxation dy- 
namics is dominated by the slow thermally activated crossing of potential 
energy barriers [3,4]. If the temperature is further reduced the relaxation time 
eventually becomes of the same order of the observation time and the sys- 
tem falls out of equilibrium since there is not enough time to cross barriers 
and equilibrate. This define the "experimental" glass transition temperature 
T g . The regime between Tmct and T g cannot be described by the MCT since 
it neglects activated processes responsible for barrier crossing. In MCT the 
relaxation time diverges at T MC t, leading to T g = Tmct-, and the dynamics 
remains confined into a single basin forever. 

The essential features of MCT for glass-forming systems are also common to 
the high temperature phase of some fully connected spin glass models [5], the 
most well known being the spherical p-spin spin glass model [6,7]. We shall call 
these models mean-field p-spin-like glass models. As a consequence at the criti- 
cal temperature T MC t, called T D in p-spin language, an ergodic to non-ergodic 
transition takes place. Below this temperature the system is dynamically con- 
fined to a metastable state (a basin) [8] since relaxation to true equilibrium 
can only take place via activated processes, absent in mean-field models. For 
these systems, nevertheless, it is known that the true equilibrium transition to 
a low temperature phase occurs below Tjj at the static critical temperature T c , 
also denoted by Ti rs b [6]. This is the analogous of the Kauzmann temperature 
Tfc for liquids. The glass transition temperature T g of real systems sits some- 
where in between T c and T D . This transition, obviously, cannot be reached 
even on infinite time in mean-field models. 

Despite these difficulties mean-field models, having the clear advantage of be- 
ing analytically tractable, have been largely used to study the properties of 
fragile glassy systems, especially between the dynamical temperature To and 
the static temperature T c . The picture that emerges is however not complete 
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since activated process cannot be captured by mean-field models. Therefore 
the relevance of mean-field results for real systems cannot be considered com- 
pletely stated. 

Only recently activated processes in mean-field-like models have been invesi- 
gated in extended numerical investigation of finite-size fully-connected p-spin- 
like models [4,9]. Comparing the results with the observed behavior of super- 
cooled liquids near T MC t we can conclude that, once activated process are 
allowed, mean-field p-spin-like models are highly valuable for a deep under- 
standing of the glass transition in real systems. 

We report the main results obtained for the Ising-spin Random Orthogonal 
Model (ROM) [10,11], defined by the Hamiltonian [10,11], H = -2 J2ij Jij °i °j 
where Oi = ±1 are N Ising spin variables, and is a iVx iV random symmetric 
orthogonal matrix with Ju = 0. For N — > oo this model has the same thermo- 
dynamic properties of the p-spin model: a dynamical transition at Tp = 0.536, 
with threshold energy per spin e t h = E t h/N = —1.87, and a static transition 
at T c = 0.25, with critical energy per spin e\ rs b = —1.936 [10,11]. 



2 Thermodynamics of Inherent Structures: How to evaluate the 
configurational entropy 

The free energy analysis (TAP) [8,11] reveals that the phase space is com- 
posed by an exponentially large (in N) number of different basins, separated 
by infinitely large (for N — > oo) barriers. Each basin is unambiguously la- 
belled by the value of the energy density e of the local minimum contained 
within it, i.e. the IS of the system. In this picture the dynamical transition is 
associated with IS having the largest basin of attraction for iV — > oo, while 
the static transition with IS with the lowest accessible free energy (vanishing 
configurational entropy) [12,8]. 

In the mean-field limit, the allowed values of e are between ei rs b and e t h- 
Solutions with e larger than e t h are unstable (saddles), while solutions with e 
smaller than e lrsb have negligible statistical weight. Moreover in the N — > oo 
limit IS with e = e t h attract most (exponentially in N) of the states and 
dominate the behavior of the system. Other IS are irrelevant for N — > oo. For 
finite N the scenario is different since not only the basins of IS with e < e t h 
acquire statistical weight, but it may happen that solutions with e > e t h and 
few negative directions (saddles with few downhill directions) become stable, 
simply because there are not enough degrees of freedom to hit them. 

To get more insight the IS-structure of finite systems we follow Stillinger and 
Weber [13] and decompose the partition sum into a sum over basins of different 
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IS and a sum within each basin. Collecting all IS with the same energy e, 
denoting with exp[Ns c (E)} de the number of IS with energy between e and 
e + de, and shifting the energy of each basin with that of the associated IS, 
the partition sum can be rewritten as [13] 

Z N (T) ~ J deexp N [-(3e + s c (e) - (3f((3, e)] (1) 

where e) can be seen as the free energy density of the system when con- 
fined in one of the basin associated with IS of energy e. The function s c (e) 
is the configurational entropy density also called complexity. From the par- 
tition function we can compute the average internal energy density u(T) = 
(e + d{(5f)/d(5) = (e(T)) + (Ae(T)). The first term is the average energy of 
the IS relevant for the thermodynamics at temperature T, while the second 
is the contribution from fluctuations inside the associated basins. In the limit 
iV — > oo only IS with e = e t h contribute and lim7v_ >00 (e(T)) = e t h for any 
T > To- For finite N, and T not too close to Tp, the thermodynamics is 
dominated by IS with e > e t h and (e(T)) > e th [4]. This is more evident from 
the (equilibrium) probability distribution of e since it is centered about (e(T)) 
indicating that IS with e ~ (e(T)) have the largest basins. This scenario has 
been also observed in real glass-forming systems [14-1 8]. 

From the knowledge of IS-energy distribution we can reconstruct the complex- 
ity s c (e) since from eq. (1) the probability that an equilibrium configuration 
at temperature T = 1/(3 lies in a basin associated with IS of energy between 
e and e + de is: P/v(e, T) = exp iV [—/3e + s c (e) — {3f(f3, e)] /Z^{T). In the 
temperature range where this applies, the curves In P N (e,T) + f3e are equal, 
except for a temperature dependent factor \nZ N (T), to s c (e) — f3f((3 ) e). If 
the e-dependence of f((3, e) can be neglected, then it is possible to superim- 
pose the curves for different temperatures, see Fig. 1 (a). The data collapse is 
rather good for e < —1.8. Above the curves cannot be superimposed anymore 
indicating that the e-dependence of e) cannot be neglected. In liquid this 
is called the anharmonic threshold [19,20]. 

Direct consequence of f(/3,e) ~ f(P) for e < —1.8 is that in this range the 
partition function can be written as the product of an intra-basin contribution 
[exp(— Nf3f)] and of a configurational contribution which depends only on the 
IS energy densities distribution. The system can then be considered as com- 
posed by two independent subsystems: the intra-basin subsystem describing 
the equilibrium when confined within basins, and the IS subsystem describing 
equilibrium via activated processes between different basins. As the temper- 
ature is lowered and/or iV increased the two processes get more separated in 
time and the separation becomes more and more accurate. A scenario typical 
of super-cooled liquids near the MCT transition [21,3]. 
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Fig. 1. (a) Configurational entropy as a function of energy. The data are from system 
sizes N = 48 (empty circle) and N = 300 (filled circle), and temperatures T = 0.4, 
0.5, 0.6, 0.7, 0.8, 0.9 and 1.0. For each curve the unknown constant has been fixed to 
maximize the overlap between the data and the theoretical result [11]. The line is the 
quadratic best-fit. (b) Configurational entropy density as a function of temperature. 
The line is the result from the best-fit of s c (e) while the symbols are the results from 
the temperature integration of d s c (T) / d (e(T)) = T _1 for N = 48 (empty circle), 
N = 300 (empty triangle) and N = 1000 (filled circle). 

3 Non-equilibrium behavior: the role of activated processes 



More informations on the IS structure can be obtained from non-equilibrium 
relaxation processes. To study the non-equilibrium dynamics we quench at 
time zero the system from an initial equilibrium configuration at temperature 
Tj > T g to a final temperature Tf < T g and study the evolution of the average 
IS energy per spin (e(t)) as function of time, see Fig.2 (b). Two different 
relaxation processes are clear seen. A first regime independent of Tf, and 
a second regime independent of both Tj and Tf. The final temperature Tf 
controls the cross-over between the two regimes. A similar behavior has been 
observed in molecular dynamics simulations of super-cooled liquids [18]. The 
two regimes are associated with different relaxation processes. In the first part 
the system has enough energy and relaxation is mainly due to path search out 
of basins through saddles of energy lower than k B Tf. This part depends only 
on the initial equilibrium temperature T, since it sets the initial phase space 
region. Different Tj leads to different power law. In particular relaxation must 
slow down as Tj decreases since we expect that lower states are surrounded by 
higher barriers, in agreement with numerical data [9]. 
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Fig. 2. Average inherent structure energy in equilibrium as a function of temperature 
(a) and as a function of time during the non-equilibrium process (6). The system 
size is N = 300, T; = 3.0 and (top to bottom) T f = 0.1, 0.2, 0.3 and 0.4 (panel b). 

During this process the system explores deeper and deeper valleys (basins) 
while decreasing its energy. The process stops when all barrier heights become 
of 0(/c#Ty). From now on the relaxation proceeds only via activated process. 
A first consequence is that lower the final temperature Tf shorter the first 
relaxation, in agreement with our findings [See figures 2]. 

The analysis of the distance between the instantaneous system state and the 
corresponding IS, counting the number of single spin flip needed to reach the 
IS, reveals that for all times the systems stays in configurations few spin flips 
away from an IS. A similar study starting from equilibrium configurations at 
temperature T e ((e(t))) evaluated comparing panels (a) and (b) of figure 2 [18] 
leads to similar numbers. We then conclude that during relaxation the aging 
system explores the same type of minima (and basins) visited in equilibrium 
at temperature T e . Direct consequence is that once the system has reached 
the activated regime there cannot be memory of the initial Tj, and all curves 
with different Tj but same Tf should collapse for large time [9]. 



4 Conclusions 

To summarize, we have shown that finite-size mean-field p-spin-like models 
are good candidates for studying the glass transition. The key point is that 
near the glass transition the thermodynamics of the systems is dominated 
by the IS distributions, therefore all systems with similar IS distributions 
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should have similar behavior. Finite-size mean-field p-spin-like models have 
the double advantage of being analytically tractable for N — > oo and easily 
simulated numerically for finite N, offering good models to analyze the glass 
transition. 
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